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HYDRODYNAMIC  TARGET  RESPONSE  TO  AN  INDUCED  SPATIAL 
INCOHERENCE-SMOOTHED  LASER  BEAM 

I.  INTRODUCTION 

One  of  the  critical  elements  for  high  gidn  laser  target  designs  is  the  high  degree  of 
symmetry  that  must  be  maintained  in  the  implosion  process.  Asymmetries  in  the  abla¬ 
tion  pressure  must  be  smaller  than  a  few  percent  if  high  gain  is  to  be  achieved.^  This 
places  a  severe  requirement  on  the  tiniformity  of  the  laser  beam  since  thermal  diffusion 
is  only  effective  at  reducing  those  spatial  nonuniformities  with  scalelengths  smaller  than 
the  separation  distance  between  the  absorption  and  the  ablation  region.^  With  short  laser 
wavelengths,  required  for  good  laser-target  coupling^,  thermal  diffusion  is  ineffective  at 
reducing  all  but  the  highest  spatial  frequency  nonuniformities. 

Recently,  three  techniques  have  been  developed  which  show  some  promise  for  reducing 
ablation  presstire  nonuniformities  down  to  fts  1%.  These  methods  are  based  on  the  con¬ 
cept  of  slicing  the  laser  light  into  small  beamlets  that  are  subsequently  integrated  at  the 
target  plane.  The  random  phase  plate  (RPP)  technique  *  slices  the  laser  light  into  many 
beamlets  with  a  transparent  phase  mask  which  randomly  imposes  a  fixed  phase  shift  of 
either  0  or  tt  on  each  of  the  beamlets.  These  beamlets  are  overlapped  onto  the  target  plane 
by  a  focusing  lens,  resulting  in  a  complicated  interference  pattern.  The  random  phase 
relationship  between  the  beamlets  remains  fixed  in  time;  as  a  result,  the  focal  interference 
pattern  remains  stationary.  In  the  smoothing  by  spectral  dispersion  (SSD)  technique,® 
the  transparent  phase  mask  is  illuminated  by  periodic  broadband  light  whose  different 
spectral  components  sweep  across  the  beam  in  cycles.  This  “color  cycling”  suppresses  the 
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interference  among  most  of  the  beamlets,  although  it  does  allow  residual  long  scalelength 
structure  due  to  interference  among  neighboring  beamlets. 


In  the  induced  spatial  incoherence  (ISI)  concept  a  brozulband  laser  beam  (bandwith 
Ai/  l/t|>«{«e)  is  sliced  into  an  array  of  small  beamlets  by  an  orthogonal  pair  of  echelon 
structures  which  impose  a  time  delay,  At,  at  each  step.  If  the  delay  increments  are  some¬ 
what  longer  than  the  optical  coherence  time,  t^  =  1/At/,  the  beamlets  become  mutually 
incoherent.  The  beamlets  are  then  overlapped  onto  the  target  by  a  lens  of  focal  length  f 
giving  a  smoothly  varying  far-field  diffraction  profile  of  total  width  2/A/jDi,  where  Di  is 
the  echelon  step  size  and  A  is  the  laser  wavelength.  At  any  instant,  the  overlapped  beam- 
lets  produce  a  complicated  interference  pattern;  however,  on  time  scales  long  compared  to 
tc,  the  interference  pattern  averages  out,  leaving  only  the  smooth  diffraction  profile.  If  the 
hydrodynamic  response  time  of  the  target  (f^)  is  much  longer  than  the  coherence  time  of 
the  laser  (tc  ^  lps)>  then  the  target  will  respond  hydrodynamically  only  to  the  smooth 
time-averaged  diffraction  profile. 

Ntimerical  simulations  of  the  time  averaged  intensity  and  ablation  pressme  profiles 
for  a  planar  target  indicate  that  an  ISI-smoothed  laser  beam  cm  restilt  in  a  peak-to-valley 
deviation  from  the  ensemble-averaged  ablation  pressure  of  less  than  2%.®  These  results  were 
obtained  with  a  3D  numerical  code  (CHAOS),  which  evaluates  the  diffraction  integral  for 
an  incident  laser  beam  of  arbitrary  aberration.  Quasi-steady  state  conditions  were  assumed 
for  the  target  and  plasma  and  thermal  smoothing  effects  were  modeled  with  a  heuristic 
“cloudy  day”,  exp(-kd),  approximation  with  an  assumed  absorption  to  ablation  distance 
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d  =  50  fim  with  0.527  fim  laser  light.*  This  relatively  good  press\ire  uniformity  would  be 
expected  as  the  dominant  mode  for  the  time  averaged  ISI-spectnun  for  the  chosen  beam 
parameters  (averaged  over  100  te)  has  a  scalelength  of  0(50  fim).  Thus  thermal  diffusion 
would  smooth  out  residual  asymmetries  with  scalelengths  of  this  order  and  smaller. 

Transverse  thermal  conduction  is  effective  at  minimiziTig  irradiation  nonuniformities 
as  long  as  there  is  sufficient  distance,  “d**,  between  the  laser  absorption  region  and  the 
ablation  re^on.  This  is  difficult  to  achieve  during  the  initial  phase  of  a  reactor-like  laser 
pulse  for  which  the  laser  intensity  is  <  lO^^W/cm^.  The  low  intensity  of  this  initial  phase 
is  reqiiired  in  order  to  keep  the  fuel  on  a  low  adiabat.  Thermal  smoothing  is  not  effective 
for  the  first  several  nanoseconds  of  the  pulse  as  the  absorption-to-ablation  distance  is 
only  a  few  microns,  or  less,  with  1/4  fxm  laser  light.  This  implies  that  the  shock  structure 
generated  during  the  start-up  phase  of  the  laser  pulse  will  mirror  the  laser  nonuniformities. 
The  shocks  produced  dxiring  the  start-up  phase  can  be  the  determining  factor  in  the  success 
or  failure  of  the  implosion  process.  The  start-up  problem  is  the  subject  which  we  address 
in  this  paper.  We  present  the  results  from  a  series  of  numerical  simulations  using  the 
FAST2D  Laser  Matter  Interaction  Model;  and,  in  particular,  we  investigate  the  impact  of 
a  shaped,  reactor-like,  ISI-smoothed  laser  beam  on  solid  DT  targets.  In  the  next  section 
we  discuss  the  numerical  model,  the  results  are  presented  in  Section  III,  followed  by  the 
summary  and  conclusions. 


3 


II.  MODEL 


The  numerical  results  are  generated  with  the  FAST2D  Laser  Matter  Interaction 
Model.  This  is  a  fully  two-dimensional,  Cartesian,  fully  compressible  hydrodynamics  code 
with  a  sliding  Eulerian  grid  with  variable  grid  spacing.  The  grid  spacing  is  typically  0.25 
fim  for  15  -  20  zones  on  either  side  of  the  ablation  layer  and  increases  uniformly  to  2  -  10 
fim  for  most  of  the  rest  of  the  gnd.  The  grids  at  the  outer  boundaries  of  the  computa¬ 
tional  mesh  are  stretched  to  200  -  300  /im.  The  grid  slides  with  the  average  velocity  of 
the  accelerated  foil  so  that  the  ablation  front  remcuns  approximately  centered  in  the  finely 
zoned  region  thus  reducing  convective  errors. 

The  time-dependent  equations  for  the  continuity  of  mass,  momentum  and  total  energy 
axe  integrated  nxunerically  with  sliding-zone  flux-corrected-  transport  (FCT).’^  FCT  solves 
the  continuity  equations  for  the  fluxes  relative  to  moving  zone  boundaries.  The  overall 
scheme  is  second-order  acctirate  and  no  artificial  viscosity  is  explicitly  required  to  maintain 
ntimerical  stability.  Mass,  momentum  and  energy  are  both  locally  and  globally  conserved. 

Thermal  conduction  is  incorporated  in  a  single-temperat\ire  model  using  a  semi- 
implicit  solution  of  the  two-dimensional,  classical,  non-linear  thermal  conduction  equation 
based  on  the  model  of  Spitzer  and  Harm.*  The  laser  energy  is  absorbed  by  classical  inverse 
bremsstrahlung  absorption  with  a  resonant  diunp  at  the  critical  density.  For  most  of  the 
results  presented  here  the  laser  light  is  at  normal  incidence.  In  order  to  investigate  the  role 
of  the  angle  of  incidence  of  the  laser  beam,  we  have  incorporated  a  ray  trace  model®  based 


4 


on  one  developed  by  Ddettrez  and  Grandjouan.^°  The  fluid  equations  are  closed  with  a 
real  equation  of  state. 

The  equations  of  state  express  the  pressure  P  =  F(e,/9)  and  the  temperature  T  = 
G(e,/)),  where  e  is  the  energy  and  p  is  the  density,  and  are  based  on  the  solid  material 
models  used  in  the  analytic  model  (ANEOS)  in  the  CHARTD  code.^^  The  model  includes 
the  effects  of  solid  material,  Fermi-temperature  and  fractional  ionization.  This  includes  the 
effect  of  a  Debye-Gruneisen  model  near  solid  density  and  the  scheme  transitions  smoothly 
from  a  Thomas-Fermi-Dirac  model  for  degenerate  electrons  at  high  density  and  low  tem¬ 
perature  to  a  Saha  equiblibrium  model  at  high  temperature  and  low  density. 

In  FAST2D,  the  ISI  laser  intensity  is  ^ven  by* 


i(v.«)  =  Wy))i’WlS(y,«)l’,  (1) 

where  (  I(y)  )  is  the  average  spatial  profile  (assumed  constant),  P(t)  is  the  slow  time 
varying  pulse  envelope  and  S(y,t)  describes  the  statistical  behavior  of  the  rapidly  varying 
stochastic  optical  field  : 

s(y.<)  =  ® 

Here,  N,  is  the  number  of  echelon  steps,  Ati  is  the  delay  increment  between  two 
adjacent  steps,  is  an  array  of  N,  statistically  independent  random  phases,  doo  is  the 
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zero  -  to  -  zero  distance  (spot  size)  and  F(t-nAti)  is  the  stochastic  complex  optical  field 


amplitude: 

F(t-nAt,)  s  ^[F,  +  iF.), 

{  F  )  =  0,  (  1  F  |>  )  =  1. 

Expanding  the  sum  over  the  number  of  echelon  steps,  we  obtain  : 


/(»,«)  =  ^  lF(t  -  nA<.)|’+ 

•  n=l 


+{FijFrk  -  FrjFik)sin[4n  -  +  (fc  -  j)^ny/doo]}- 


Fr  and  Fj  are  independent  arrays  of  real  munbers  satisfying  Gaussian  statistics,  and 

(  Fr  )  =  (  F.  )  =  0;  (Fr)rm.  =  (F.)rm.  =  1. 

These  arrays  are  regenerated  every  numerical  time  step  and  the  numerical  <.ime  step 
is  forced  to  remain  <  tc.  As  noted  above,  the  spot  size  is  given  by  doo  =  2  f  A/Di,  where 
f  is  the  focal  length,  A  is  the  laser  wavelength  and  Di  is  the  step  size.  The  length  of 
the  n**  transverse  mode  in  the  ISI  spectrum  is  given  by  An  =  f  A/(nDi).  The  longest 
mode  is  doo/2  and  the  shortest  mode  is  doo/(2N,).  FAST2D  employs  periodic  boundary 
conditions  in  the  direction  transverse  to  the  laser  propagation,  thus,  computationally,  the 
longest  mode  represented  is  doo/2;  i.e.,  half  of  the  spot  size.  The  laser  beam  is  assumed 
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to  have  a  flat  spatial  proflle;  i.e.,  (/(y)}  =  constant.  The  numerical  results  are  presented 
in  the  next  section. 

III.  RESULTS 

For  most  of  the  results  presented  here,  we  modeled  the  laser  proflle  with  a  20  x  20 
echelon  array  with  doo  ~  320  pm  and  a  laser  coherence  time  te  =  1  ps.  The  laser  wavelength 
is  0.264  pm.  This  corresponds  to  a  lens  with  an  Fnvmter  =  32.  The  laser  spot  is  flat  over 
the  distance  doo-  The  transverse  dimension  is  160  pm,  represented  by  80  zones,  and  there 
are  120  zones  in  the  longitudinal  direction.  The  ISI  mode  spectrum  ranges  between  8 
pm  and  160  pm.  This  compares  favorably  with  the  envisioned  reactor-like  parameters  ; 
Fnum&er  »  30,  doo  =  600  pm  and  a  60  x  60  echelon  array.  These  conditions  produce  an 
ISI  spectrum  ranging  between  5  pm  and  300  pm. 

The  instantaneous  incident  ISI  laser  profile  exhibits  large  amplitude  randomly  fluc¬ 
tuating  structures  as  shown  in  Figure  1.  Here  is  plotted  the  instantaneous  incident  laser 
intensity  (in  arbitrary  units)  versus  the  transverse  dimension  (in  microns)  at  four  arbitrary 
times.  The  long  time  average  behavior  (100  te)  still  exhibits  some  residual  structure  as 
shown  in  Figure  2.  After  100  ps,  the  laser  intensity  has  a  root-mean-square  deviation 
<^RMS  =  0.10  with  a  maxmimum  amplitude  variation  of  1.7:1,  or  ±25%.  The  RMS  devi¬ 
ation  for  the  incident  light  follows  a  (tc/<)‘^*  dependence.®  This  is  illustrated  in  Figure  3 
which  shows  the  time  averaged  incident  laser  profile  after  500  tc  and  after  1000  tc-  After 
1  ns,  a  RMS  =  0.037  with  a  maximtxm  amplitude  variation  of  ±8.6% 
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These  residual  structures  can  be  reduced  in  magnitude  by  lateral  thermal  smoothing 
as  long  as  the  distance  between  the  ablation  region  and  the  laser  absorption  region  is 
larger  than  the  scalelengths  of  the  residual  laser  asymmetries.  To  test  this  hypothesis, 
we  generated  a  moderately  long  scalelength  plasma  with  a  “perfect”  laser  beam  before 
impacting  the  target  with  an  ISI-smoothed  laser  beam.  The  laser  pulse  has  a  2  ns  Gaussian 
rise  to  5  X  10^^  W/cm*  after  which  it  has  a  power  law  dependence  (t“®/^)  which  peaks  at  5 
X  10^^  W/cm^  at  13  ns  and  is  then  held  constant.  The  laser  light  is  spatially  uniform  imtil 
6  ns,  at  which  time  the  ISI  spectnim  is  “turned  on”.  At  this  time  the  distance  between 
the  ablation  and  critical  surfaces  is  fs  1/im;  the  distance  between  the  ablation  and  1/4 
critical  surfaures  is  15^m.  Nearly  all  of  the  laser  energy  is  absorbed  between  the  1/4 
critical  density  and  critical  density  surfaces.  Figure  4  shows  the  isodensity  conto\irs  and  a 
2D  perspective  plot  of  the  mass  density  of  an  80  fim  thick  CH  target  at  16  ns.  The  target 
is  quite  uniform  well  into  the  rapid  acceleration  phase.  The  mass  perturbation  at  16  ns 
is  Am/m  »  0.2%.  It  is  evident  that  transverse  thermal  conduction  effectively  reduces  the 
long  time  average  residual  structure  inherent  in  the  ISI  l2iser  beam  profile  when  there  is 
suflBcient  separation  between  the  laser  absorption  region  and  the  ablation  region. 

These  density  profiles  are  to  be  contrasted  with  the  case  of  a  “£rozen-in”  4:1  sinusoidal 
laser  asynunetry  which  is  illustrated  in  Figure  5.  The  initial  conditions  axe  the  same  as  the 
ISI  case;  i.e.,  the  laser  profile  is  uniform  until  6  ns  at  which  time  the  target  is  illuminated 
with  the  asymmetric  beam.  At  9.5  ns,  3.5  ns  after  the  asymmetric  beam  is  “turned  on  ,  the 
target  is  severely  distorted  with  a  13  nm  depression  which  mirrors  the  laser  asymmetry. 
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Note  that  this  is  solely  due  to  the  differential  mass  ablation  rate  and  ablation  pressvire 
asymmetry.  There  is  very  little  target  acceleration  at  this  point  in  the  laser  pulse  and  the 
Rayleigh-Taylor^^’^^  (RT)  instability  is  not  playing  a  role  here. 

A  similar  result  is  obtained  with  an  ISI-like  beam  with  zero  bandwidth.  In  this  case, 
at  6  ns,  the  uniform  laser  profile  shifts  to  a  “frozen-in”  ISI-profile.  A  plot  of  the  incident 
intensity  is  shown  in  Fig\ire  6.  The  target  density  contours  are  shown  in  Figure  7  at  11  ns. 
The  target  is  markedly  distorted  at  this  time  and  a  portion  of  this  growth,  »  1  e-folding, 
is  due  to  the  RT  instability  as  the  target  has  begim  to  accelerate  at  this  point  in  the  pulse. 
At  6  ns,  the  ablation-to-absorption  distance  is  not  large  enough  to  thermally  smooth  out 
the  large  scalelength  asymmetries  resulting  from  the  zero  bandwidth  ISI  case. 

The  ability  of  the  ISI  concept  to  result  in  a  smooth  ablation  pressure  is  based  on  a 
combination  of  temporal  smoothing  (over  lOO’s  of  tc)  and  transverse  thermal  smoothing. 
Transverse  thermal  conduction  is  effective  at  minimizing  the  long  time  average  residual 
structure  inherent  in  the  ISI-smoothed  laser  beam  as  long  as  the  sotmd  transit  time  between 
the  laser  absorption  repon  and  the  ablation  front  is  greater  than  the  ISI  averaging  time 
(100  -  200  tc).  Unfortxmately,  this  necessitates  a  long  scalelength,  uniform  plasma  in  the 
laser  absorption  region.  This  is  difficult  to  achieve  with  a  reactor-like  shaped  laser  pulse 
which  may  entail  a  long,  low  intensity  “foot”  prior  to  the  main  driving  beam.  The  laser 
intensity  can  be  <  10^®W/cm*  for  nearly  80%  of  the  temporal  length  of  the  pulse.  A 
canonical  power  law  shaped  laser  pulse  is  shown  in  Figure  8.  After  a  1/2  ns  Gaussian  rise 
to  10^*  W/cm*,  the  pulse  has  a  t“*/^  dependence  to  3  x  10^*  W/cm^  at  15  ns  after  which 


the  laser  intensity  is  held  fixed.  This  ptilse  is  designed  so  that  the  first  shock  breaks  out 
through  the  rear  of  the  target  as  the  laser  intensity  reaches  its  maximum  value. 

The  severity  of  the  thermal  smoothing  problem  is  illustrated  in  Figure  9,  where  the 
distances  between  the  ablation  surface  and  the  critical  surface  (solid  line)  and  the  ablation 
surface  and  the  quarter  critical  siirface  (dashed  line)  are  plotted  as  a  function  of  time  for 
the  above  laser  pulse.  For  2/3  of  the  laser  pulse,  the  critical  surface  is  less  than  1/2  /im 
from  the  target  surface. 

The  initial  pulse  sends  an  w  2/3  Mb  shock  through  the  330  fim  thick  DT  target  which 
keeps  the  fuel  on  a  low  adiabat.  The  ablation  pressure  history  and  the  peak  target  dexisity 
history  are  shown  in  Figure  10.  The  target  density  is  initially  compressed  by  a  factor  of 
4,  from  0.2  to  nearly  0.8  gm/cm^.  For  a  strong  shock,  the  compression  ratio  is  given  by 
(7  +  l)/(7  -  1)  where  7  is  the  ratio  of  specific  heats,  7  «  5/3  for  DT.  At  the  peak  of  the 
laser  pulse,  the  target  density  has  been  compressed  by  a  factor  of  over  30.  The  ablation 
pressure  at  this  time  is  45  Mbar.  The  velocity  history  of  the  target  is  shown  in  Figure  11. 
The  initial  shock  imparts  a  velocity  of  1.5  x  10*  cm/s  to  the  target.  The  target  velocity 
remains  nearly  constant  imtil  12  -  13  ns,  at  which  time  the  target  begins  to  rapidly 
accelerate  as  a  result  of  the  increasing  laser  intensity.  The  target  velocity  is  7.5  x  10*  cm/s 
at  15  ns  and  the  acceleration  is  8.5  x  10'*  cm/s*.  At  this  point  the  target  has  traveled 
300  /im.  It  is  the  mass  variations  near  the  peak  of  the  laser  pulse  which  will  provide  the 
seeds  for  the  RT  growth  during  the  acceleration  phase.  For  a  target  to  implode  uniformly, 
the  mass  variations  at  this  point  must  be  ^  1%.  Note  that  the  target  attains  a  sizeable 
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acceleration  between  the  12  -  15  ns  time  interval  and  any  residual  mass  perturbations  can 
undergo  1-2  e-foldings  of  RT  growth  during  this  time  period.  For  the  results  presented 
below,  the  ISI  spectrum  is  in  effect  from  the  beginning  of  the  pulse. 

The  isodensity  contours  for  a  330  fim  thick  DT  target  illuminated  with  a  0.264  fim 
ISI  laser  beam  is  shown  in  Figure  12  at  t  =  12  ns.  At  this  time  the  target  mass  variation 
{A(pr)/{pr})  is  22%.  Early  in  the  laser  pulse,  the  laser  energy  is  deposited  directly  onto 
the  target  surface.  The  highly  nonuniform  nature  of  the  ISI  beam  results  in  spatially 
and  temporally  nonuniform  shocks  traversing  the  target.  The  target  surface,  itself,  is 
perturbed  and  as  the  shock  passes  through  the  target  it  leaves  residual  mass  perturbations 
in  its  wake.  This  is  apparent  from  Figure  13  where  plots  of  the  target  mass  variation 
(A(pr)/{pr))  at  0.4,  0.6,  0.8  and  1.0  ns  show  rapid  variations  in  the  mass  perturbations. 
The  same  phenomenon  is  discemable  in  2D  perspective  plots  of  the  mass  density.  Figure 
14.  As  the  shock  traverses  the  dense  target,  pertmbations  on  the  shock  front  eventually 
dissipate  and  the  shock  front  straightens  out.  This  occurs  because,  relative  to  the  average 
velocity  of  the  shock  front,  the  leading  portion  of  the  shock,  the  crest,  is  divergent  and  thus 
slows  down;  whereas,  the  trailing  portion,  the  trough,  is  convergent  and  thus  speeds  up. 
In  addition,  once  the  sound  transit  time  between  the  absorption  region  and  the  ablation 
region  becomes  greater  than  the  ISI  averaging  time,  the  rapidly  varying  ISI  laser  beam  no 
longer  pertxirbs  the  target.  The  residual  mass  variations,  however,  persist.  Mass  density 
plots  at  3.0  and  3.2  ns  shown  little  difference  in  the  mass  variation  on  the  target  smface, 
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see  Figure  15.  These  residual  variations  grow  linearly  in  time,  in  a  Richtmyer-Meshkov^^- 
like  manner,  as  the  target  drifts  with  a  nearly  constant  velocity.  Figure  16  shows  a  2D 
surfcu;e  plot  of  the  integrated  mass  density  (pr)  as  a  function  of  time,  pr  is  integrated  from 
the  rear  of  the  target  to  the  ablation  front  for  each  transverse  coordinate.  The  residual 
perturbations  begin  to  grow  after  about  1  ns.  Apr  between  line-outs  taken  along  paths  1 
and  2  and  along  paths  2  and  3  (see  arrows  in  Figure  16)  are  plotted  as  a  function  of  time 
in  Figure  17.  The  initial  growth  is  linear  in  time. 

That  this  should  be  the  case  is  clear  from  the  theoretical  expression  for  RM  growth*®. 


T}  =  i2irlX)rf^AViil  -  Ui/U.),  (4) 

where  rjo  is  the  initial  amplitude,  A  is  the  wavelength,  U,  is  the  velocity  of  the  shock  wave, 
Ui  is  the  velocity  imparted  to  the  interface,  the  Atwood  number  A’  =  (p!j.  —  pL)/(pV  +P-)> 
where  p+,pL  are  the  densities  at  the  interface  just  after  passage  of  the  shock.  For  the 
parameters  investigated  here,  A'  is  0(1),  and 


AT}/r}o  =  {2n/X)Ui{l  -  UilU,)iXt.  (5) 

The  appropriate  time  increment  is  the  shock  breakout  time,  At  =  Ax/C/^*,  where  Ai 
is  the  target  thickness.  From  the  mass  conservation  equation,  Ui/Ut  =  (1—  p<i/p«)» 
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where  pd{Pu)  is  the  density  downstream  (upstream)  of  the  shock.  For  a  strong  shock, 
Pd/Pu  =  (7  -  l)/(7  +  1)  =  1/4,  for  DT,  and 

=  37r/(8A)Ai.  (6) 

Thus  for  perturbation  wavelengths  of  40  -  50  pm  and  target  thicknesses  of  300  pm,  sTna]! 
initial  perttirbations  can  grow  an  order  of  magnitude  in  size  solely  as  a  result  of  shock  driven 
growth,  well  before  the  target  acceleration  phase.  This  is  eqmvalent  to  over  2  e-foldings  of 
growth.  This  growth  will  occur  whether  the  perturbations  stem  from  laser  imperfections 
or  target  surface  imperfections^®.  The  problem  is  exacerbated  by  thick  targets  and  strong 
shocks. 

Further  evidence  for  the  claim  that  the  initial  shock  structure  is  the  determining 
factor  in  the  subsequent  evolution  of  the  target  can  be  gained  from  the  case  in  which 
the  ISI  spectrum  is  “turned  off”  after  3  ns;  i.e.,  the  laser  profile  is  modeled  as  being 
perfectly  uniform  for  t  >  3ns.  There  is  essentially  no  difference  in  the  amplitude  of  the 
mass  perturbations,  later  in  time,  for  this  case  as  compared  to  the  case  in  whidi  the  ISI 
spectrum  is  “on”  throughout  the  whole  pulse.  In  addition,  a  case  was  nm  in  which  the  ISI 
spectrum  is  averaged  over  500  tc  and  then  held  fixed  (see  Figxire  3a)  during  the  shaped 
power  law  laser  pulse.  After  3  ns,  the  amplitude  of  the  mass  pertmbation  is  <  10%  smaller 
than  for  the  full  ISI  spectrum  case. 
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Although  it  appears  that  a  long  pulse  ISI  laser  beam  may  severely  pertmb  an  ICF 
target,  this  should  be  contrasted  with  the  impact  of  a  temporally  equivalent  laser  pulse 
with  a  ‘%ozen-in”  4:1  sinusoidal  asymmetry.  In  this  case,  the  laser  as}nnmetry  is  imposed 
at  t  =  0  and  the  mass  variation  is  greater  than  2:1  at  8  ns. 

Several  target  designs  were  investigated  in  an  att'-.tupt  to  mitigate  the  impact  of  the 
initial  shock  structure  stemming  from  the  early  temporal  phase  of  the  ISI  laser  beam. 
In  one  case,  a  180  /im  thick  solid  DT  target  with  a  100  ^lm  thick  ‘Yoam-like”  DT  layer 
was  vised.  This  “foam-like”  layer  does  not  contain  any  voids  and  pore  comp2u:tion  is  not 
modeled.  The  “foam-hke”  layer  is  simply  low  density  DT,  p/oam  =  O-OSgm/cm^.  This 
layer  reduces  the  perturbation  on  the  shock  front  to  some  degree  before  the  shock  impacts 
the  frozen  DT  target.  The  mass  variation  at  15  ns  is  21%  with  0.264  pm  ISI  laser  light. 

The  ablation  to  absorption  distance  can  be  increased  by  using  longer  wavelength  laser 
light  which  would  also  result  in  an  increase  in  the  amount  of  thermal  smoothing.  However, 
longer  wavelength  laser  light  implies  poorer  iMer-target  coupling  and  reduced  efficiency.® 
By  combining  long  wavelength  laser  light  at  the  beginning  of  the  pulse  with  short  wave¬ 
length  laser  light  during  the  compression  and  acceleration  phases,  thermal  smoothing  would 
be  enhanced  during  the  early  phase  of  the  pulse  and  laser-target  coupling  would  be  en¬ 
hanced  during  the  acceleration  phase.  Using  the  same  multilayer  target  discxissed  above, 
the  first  8  ns  of  the  laser  beam  is  modeled  as  1.054  pm  ISI  laser  light  (tc=l  ps)  which  then 
smoothly  evolves,  over  1.5  ns,  to  0.264  pm  ISI  laser  light  (tc=l  ps)  for  the  rest  of  the  laser 
pvilse.  At  15  ns,  A(pr)/(pr)  =  4%.  The  perturbation  amplitude  at  this  time  includes  « 
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1.5  e-foldings  of  RT  growth.  At  12  ns,  A(pr)/{pr)  =  0.9  %.  Modeling  the  first  8  ns  of  the 
1.054  fim  beam  as  being  perfectly  uniform  results  in  a  mass  variation  of  only  1.8%  at  15 
ns  (0.4%  at  12  ns).  These  perturbation  levels  are  still  too  large  to  enable  the  target  to 
implode  uniformly.  The  isodensity  contours  for  the  second  of  these  two  cases  is  shown  in 
Figure  18. 

The  ablation  to  absorption  distance  also  depends  on  the  angle  of  incidence  of  the  laser. 
The  position  of  the  tiiming  point  for  the  laser  light  is  given  by  rt  =  te  +  2L  /n(cos(^), 
where  rc  is  the  position  of  the  critical  surfsice,  L  is  the  density  scalelength  and  6  is  the  angle 
of  incidence  of  the  laser.  Since  it  is  the  early  time  impact  of  the  laser  which  is  perturbing 
the  target,  L  must  be  nonzero  in  order  that  there  be  any  increase  in  thermal  smoothing. 
In  the  absence  of  a  “perfect”  laser  beam,  a  uniform  blowoff  plasma  could  be  generated  by 
coating  the  target  with  a  thin  (0(few  100  A)  high-Z  material  and  striking  the  target  with 
a  very  short  burst  of  x-rays.  We  have  approximated  this  interaction  by  depositing  1/4  J  of 
energy  directly  onto  the  stirface  of  a  frozen  DT  target  during  a  250  ps  FWHM  pulse.  The 
target  is  330  /xm  thick.  This  energy  pulse  produces  a  very  weak  shock,  O(10’s  of  kbar), 
and  at  10  ns  after  the  “x-ray  fiash”,  the  density  scalelength,  at  critical  density  for  1/4  pm 
laser  light,  is  w  100pm.  The  density  profile  at  this  time  is  shown  in  Figure  19.  The  target 
is  then  illuminated  with  a  0.264  pm  ISI  -  smoothed  laser  beam  at  a  30*  angle  of  incidence 
(with  respect  to  the  normal).  For  this  case,  the  laser  pulse  has  a  5  ns  Gaussian  rise  to  an 
intensity  of  10^^  W/cm*  which  then  evolves  into  a  power  law  pulse  which  peaks  at  3  x  10*^ 
W/cm  *  at  15  ns.  A  plot  of  the  rays  incident  at  an  angle  of  30“  at  the  peak  of  the  laser 
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ptilse  is  shown  in  Figure  20.  The  turning  point  of  the  rays  is  just  outside  of  the  critical 
surface  and  the  plasma  continues  to  absorb  energy  &om  the  rays  tmtil  each  ray  contains  < 
1%  of  its  initial  energy.  The  maximum  mass  variation  of  the  target,  A(/)r)/(pr),  is  4.8  % 
at  15  ns.  The  isodensity  contours  at  15  ns  are  shown  in  Figure  21.  The  mass  perturbations 
are  larger  than  the  layered,  multiple  laser  wavelength  cases,  but  significantly  smaller  than 
the  nonlayered,  normal  incident  cases. 

IV.  SUMMARY  AND  CONCLUSIONS 

We  have  niimeiically  simulated  the  hydrodyniunic  response  of  a  series  of  thick,  frozen 
DT  targets  directly  illuminated  with  an  ISI-smoothed,  ICF  re8u:tor-like,  shaped  laser  pulse. 
Studies  of  the  ‘‘start-up”  problem  and  the  impact  of  the  first  shock  on  subsequent  target 
development  are  important  for  both  target  design  and  pulse  shaping.  The  results  indicate 
that  the  first  shock  can  be  the  determining  factor  in  the  success  or  failure  of  the  implosion 
process.  This  is  the  case  whether  the  perturbations  on  the  shock  front  stem  from  laser  im¬ 
perfections,  as  investigated  here,  or  stirface  imperfections.^*  A  pertxirbed  shock  front  leaves 
residual  mass  variations  in  its  wake  as  it  transits  the  target.  These  mass  variations  are  the 
seeds  for  the  subsequent  Rayleigh-Taylor  growth  which  will  occur  during  the  acceleration 
phase  of  the  laser  pulse.  For  a  strong  shock  and  thick  targets,  these  mass  variations  can 
undergo  «  2  e-foldings  of  growth  before  the  target  begins  the  main  acceleration  phase. 

The  ability  of  the  ISI  technique  to  provide  nearly  uniform  ablation  pressures  is  con- 
tigen!.  on  both  temporal  averaging,  over  lOO’s  of  coherence  times,  and  thermal  diffusion. 
With  a  moderate  scalelength,  uniform,  preformed  plasma,  an  ISI-smoothed  laser  beam  can 


accelerate  the  target  uniformly.  Under  similar  initial  conditions,  asymmetric  laser  beams 
and  beams  with  echelon  steps  but  zero  bandwidth  can  severely  perturb  the  target.  In  the 
latter  two  cases  there  is  no  temporal  avera^ng,  and  thermal  diffusion  is  not  large  enough 
to  smooth  out  the  intensity  perturbations. 

In  the  absence  of  a  ‘^perfect”  laser  beam  and/or  a  uniform  preformed  plasma,  some 
means  may  be  needed  to  mitigate  the  impact  of  the  first  shock  firom  a  shaped  reactor-like 
laser  pulse.  We  have  shown  that  cold,  low  density  layers,  multiple  wavelength  lasers  and 
shallow  angles  of  incidence  can  all  significantly  reduce  the  perturbations  stemming  from 
an  imperfect  laser  beam.  Some  combination  of  the  above  tsurget  designs  may  be  needed 
in  order  to  uniformly  implode  a  thick  IGF  capsule.  It  should  be  noted  that  these  designs 
may  reduce  the  overall  efficiency  of  the  implosion  process.  Finally,  it  may  be  possible  to 
compress  the  target  adiabatically  •  without  any  shodcs.  This  could  entail  a  very  long  (  50 
-  100  ns),  slowly  rising  laser  pulse.  A  laser  pulse  of  this  shape  would  generate  kilobar-like 
pressures  in  the  target,  and  thus  compressive  stress,  shear  stress,  melt  and  pore  compaction, 
in  the  case  of  foams,  become  important.  We  are  presently  modifying  our  EOS  model  in 
order  to  more  accurately  simiUate  these  physical  phenomena. 
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.  Plot  of  the  instantaneous  incident  ISI  laser  intensity  (arbitrary  units)  versus  the  trans¬ 
verse  dimension  (in  microns)  at  four  different  times. 
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2.  Plot  of  the  laser  intensity  which  has  been  averaged  over  100  coherence  times.  The  rms 

variation  is  10%. 
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5.  Isodensity  contoxirs  of  an  80  nm  thick  CH  target  impacted  with  a  laser  with  an  imposed 
4:1  sinusoidal  asymmetry.  The  time  is  9.5  ns,  the  tremsverse  dimension  is  160  fim  and  the 

horizontal  dimension  is  60  um. 
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6.  A  plot  of  the  incident  intensity  for  an  ISI  beam  with  zero  bandwidth. 
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8.  Temporal  variation  of  a  shaped  power  law  leiser  pulse. 
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9.  Plot  of  the  distance  between  the  ablation  surface  and  the  critical  surface  (solid  line)  and 
the  distance  between  the  ablation  surface  and  the  quarter  critical  surface  (dashed  line),  in 
microns,  versvis  time  for  the  shaped  power  law  laser  pvilse  illustrated  in  Figure  8. 
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12.  Isodensity  contours  at  12  ns  for  a  330  ftm  thick  DT  target  illuminated  with  a  0.264 
fim  ISI  laser  beam  with  a  power  law  time  dependence.  The  horizontal  dimension  is  124 

ftm  and  the  transverse  dimension  is  160  /xm. 
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13.  Plots  of  the  total  mass  variation  AAf/A/  versus  the  transverse  dimension  (in  microns) 
at  0.4,  0.6,  0.8, 1.0  ns  showing  the  rapid  oscillations  in  the  mass  perturbations  as  the  shock 
impacts  the  target.  AM/M  has  been  normalized  to  1  for  each  time. 


14.  2D  perspective  plots  of  the  target  mass  density  at  the  same  4  times  as  in  Figiire  13. 
The  laser  is  impacting  the  target  from  the  front.  The  vertical  dimension  is  linear. 


in  the  mass  variation  at  this  time. 
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16.  A  2D  surface  plot  of  the  integrated  mass  density  pr  versus  the  transverse  target 

dimension  (microns)  as  a  function  of  time  (ns). 
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18.  Isodensity  contours  at  15  ns  for  180  frozen  DT  target  with  a  100  fim  low  density 
layer  illuminated  with  “perfect”  1  laser  light  between  0  -  8  ns  and  ISI-smoothed  1  /4 

ftm  laser  light  between  8  -  15  ns. 


37 


X(mtcrorw) 

20.  Plot  of  the  rays  incident  at  an  angle  of  30®  at  the  peak  of  the  laser  pulse.  The  horizontal 

and  vertical  dimensions  are  in  microns. 
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21.  Isodensity  contours  for  the  30®  angle  of  incidence  case  at  the  peak  of  the  laser  pulse. 
The  horizontal  and  vertical  dimensions  are  both  160  nm. 


